甘蔗 耐寒 相关 miRNAs 的 生物 信息 学 分 析 及 其 对 基因 预测 
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摘要 : 为 了 解 不 同 基 因 型 甘蔗 (Saccparzz officinarum) 响应 低温 胁迫 的 分 子 机 制 ， 挖 掘 厂 
寒 相关 的 miRNA 及 相关 靶 基 因 ， 该 文 以 低温 胁迫 4‘C 处 理 24 h 后 的 3 个 不 同 耐寒 性 甘蔗 
品种 的 叶片 为 材料 进行 Tllumina HiSeqTM2000 高 通 量 测序 ， 构 建 18 个 低温 胁迫 前 后 sRNA 
文库 。 结 果 表 明 :(1) 共 获得 分 属于 84 个 家 族 的 已 知 miRNA 322 个 及 预测 得 到 110 个 新 miRNA， 
并 在 已 知 mi RNAs 中 筛选 出 差异 表达 miRNA100 个 (61 个 上 调 ，39 个 下 调 ) ， 新 miRNA 中 筛选 
出 差异 表达 miRNA37(15 个 上 调 ，22 个 下 调 ) (2) 利用 psRNATarget. TargetFinder 和 
Tapirhybrid 软件 对 所 获得 的 差异 表达 miRNA 进行 靶 基 因 预 测 ， 预 测 得 到 1844 个 靶 基 因 并 
进行 GO 分 析 揭 示 了 这 些 间 基因 的 3 个 主要 功能 类 别 ， 即 分 子 功能 、 细 胞 组 分 与 生物 过 程 。 
(3) 为 验证 高 通 量 测序 数据 的 可 靠 性 ， 筛 选 14 个 miRNA 及 其 靶 基 因 进 行 qRT-PCR 验证 , 22 
明 这 些 miRNA 均 被 检测 发 现 且 大 部 表达 结果 与 测序 结果 一 致 。(4) 还 鉴定 出 部 分 差异 表达 
iRNA 的 靶 基 因 ， 这 些 基 因 参 与 植物 生长 、 发 育 及 低温 胁迫 反应 。 综 上 结果 表明 ， 耐 寒 型 甘 
RETE PI miRNA 直接 或 者 间接 作用 靶 基 因 实 现 表达 调控 相关 代谢 途径 , 对 其 重要 农艺 性 状 均 起 
着 关键 的 调控 作用 。 
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Bioinformatics analysis of microRNAs and prediction of target genes associated with cold 


tolerance in sugarcane 
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Ping, Ou Kewei, Lu Yefei, Nong Zemei, Tang Huanwei, Long Shengfeng 

(Guangxi Subtropical Crops Research Institute, Nanning 530001, China) 
Abstract: To identify the molecular mechanisms of on sugarcane (Saccharum officinarum ) 
responding to cold stress, and find the miRNAs and target genes related to cold tolerance, the 
leaves of different sugarcane genotypes with different cold tolerance treated for 24 hat 4 °C were 
sampled as the materials for high-through transcriptome sequencing with IlluminaHiSeqTM 2000, 


and 18 sRNA libraries before and after cold stress were constructed. The results were as follows: 
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(1) Total 322 known miRNAs of 84 families were discovered, 110 new miRNAs were predicted. 
100 differentially expressed miRNAs were screened out from the known miRNAs (61 
up-regulated, 39 down-regulated) and 37 differentially expressed miRNAs (15 up-regulated, 22 
down-regulated) were screened out from the new miRNAs. (2) 1 844 target genes were predicted 
by using psRNATarget, TargetFinder and Tapirhybrid software. Three main functional categories 
of these target genes were revealed via the functional analysis of gene ontology, namely molecular 
function, cellular component and biological process. (3) In order to verify the high-throughput 
sequencing data, 14 miRNAs and their target genes were selected for qRT-PCR analysis, which 
showed that most of the expressions of the 14 miRNAs were detected and consistent with the 
sequencing results. (4) Some miRNA target genes were identified, which involved in plant 
growth, development, and cold stress responses. This study shows that miRNA directly or 
indirectly regulates the expression of target genes in cold-tolerant sugarcane, and plays a key role 


in regulating the important agronomic traits. 
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MicroRNAs 是 一 类 20-24 个 碱 基 的 小 分 子 内 源 性 非 编码 RNA, 
具有 高 度 的 保守 性 、 时 序 性 和 组 织 特异 性 (Ambros, 2004), miRNA = 
要 在 转录 后 水 平 上 通过 介 导 靶 基 因 mRNA 的 切割 或 抑制 翻译 来 调节 
基因 的 表达 ,在 生物 体 代谢 过 程 中 起 到 多 种 调控 作用 ,诸如 参与 调控 
植物 器 官 的 形态 建成 (Sunkar, 2004, 2012)、 生 长 发 育 (Thiebaut, et al., 
2012)、 激 素 分 泌 、 信 号 转 导 以 及 对 外 界 环境 胁迫 (Xiong, 2003) 的 响应 
过 程 。 在 拟 南 芥 和 水 稻 响 应 低温 的 研究 中 , 发 现 miR-167、 miR-169, 
miR-319 和 miR-171 等 miRNA 家 族 在 低温 响应 具有 重要 的 生物 学 作 
用 (Sunkar 2005; Wang et al., 2010)， 其 中 miR-169 的 靶 基 因为 低温 诱 
导 的 重要 基因 CBF (Sunkar et al. , 2007); 还 在 水 稻 响 应 低温 胁迫 的 研 
究 中 ,证明 miRNA-319 和 miRNA-171 的 靶 基 因 属 于 MYB 类 的 转录 
因子 ， 且 两 者 的 表达 量 互 为 消长 的 关系 ， 进 而 说 明 mRNA 在 水 稻 耐 
低温 途径 中 所 起 的 调控 作用 (Sunkar et al., 2005; Lü et al., 2010)。 因 此 ， 


准确 高 效 分 离 和 鉴定 miRNA АВА MTENDE, SETS 
地 掌握 miRNA 在 植物 搞 逆 胁迫 过 程 中 的 调控 机 制 ， 是 目前 植物 
miRNA 研究 领域 的 重要 内 容 。 

甘 蕊 是 热带 、 亚 热带 地 区 重要 的 糖 料 作物 ,低温 是 限制 其 扩大 种 
植 区 域 和 实现 高 产 稳产 的 重要 因素 之 一 ( 吴 棉 国 等 , 2010)。 出 现 的 偶 
发 性 大 学 围 冰 冻雨 雪 或 严重 的 祖 冻 天 气 将 导致 我 国 甘蔗 主 产 区 遭受 
巨大 的 经 济 损失 ( 何 燕 等 , 2009; 匡 昭 敏 等 , 2009);， 还 有 广西 地 理 环境 
特殊 ,在 冬季 持续 低温 以 及 春季 “倒春寒 ”形成 的 阴雨 霜冻 频 发 以 致 
大 面积 甘蔗 受到 冷害 ， 具 体 表现 在 甘蔗 叶片 枯萎 、 茎 杆 坏死 以 致 蔗糖 
含量 、 甘 蔗 减 产 下 降 (李杨 瑞 等 , 2011)。 有 研究 采集 甘蔗 主 产 市 县 40 
多 年 的 气象 数据 , 结合 甘蔗 种 植 面积 、 甘 莽 产 量 和 蔗糖 产量 等 指标 进 
行 分 析 发 现 低温 神 头 灾害 是 影响 甘蔗 种 植 、 甘蔗 生产 和 蔗糖 产量 的 主 
要 环境 因子 ( 何 燕 等 , 2008; 古 丽 等 , 2011D)。 因 此 ， 了 解 甘蔗 的 耐 低 温 
调控 机 制 是 培育 耐寒 性 强 或 适合 我 国 热 囊 北 缘 气 候 的 甘蔗 品种 的 前 
提 。 

本 研究 在 观察 不 同 甘 蔗 品 种 的 田间 农艺 性 状 和 模拟 低温 胁迫 的 
生理 生化 研究 基础 上 , 筛选 出 抗 寒 能 力 较 强 的 甘蔗 品种 作为 研究 材料 ， 
利用 高 通 量 测序 技术 及 生物 信息 学 方法 , 获得 与 甘蔗 响应 低温 胁迫 相 
天 的 miRNA, 分 析 miRNA 的 差异 表达 情况 , 明确 miRNAs SERA 
的 作用 关系 , 并 对 预测 所 得 的 靶 基 因 进 行 基因 本 体 (gene ontology, GO) 
分 析 , 挖掘 低温 胁迫 应 答 基因 ， 为 选 育 耐寒 性 强 的 优 展 甘蔗 新 品种 提 
供 理论 依据 和 技术 支撑 。 


1 材料 与 方法 
1.1 材料 

以 广西 农 科 院 甘蔗 研究 所 选 育 的 “ 桂 糖 28 5 (GT28). Гр 
区 主 栽 品种 “新 台糖 22 号 ”(ROC22) 和 广西 亚 热 市 作物 研究 所 新 选 
ЕН 桂 热 2 号 ”品种 (GR2) 为 实验 材料 。 
1.2 方法 
1.2.1 实验 设计 与 管理 

选择 无 病虫害 、 莽 茎 大 小 均匀 的 种 茎 切 成 单 芽 段 ， 先 用 清水 冲洗 
种 茎 并 用 干 布 擦洗 干净 ,接着 用 50% 多 菌 灵 可 湿性 粉剂 1 000 倍 液 浸 
12h 消毒 ， 再 用 蒸馏 水 浸 洗 种 茎 1 min 后 ， 用 蒸馏 水 浸 洗 拧 干 的 棉布 
将 种 茎 包 好 并 做 好 标记 ,然后 用 橡皮 筋 扎 好 放 进 温度 为 25“C 的 恒温 
箱 催 车。 当 种 茎 萌芽 并 长 出 幼 根 时 ,将 其 移植 到 装 有 和 营养 土 的 塑料 伪 
并 做 好 标记 。 每 倪 1 ЗЕ, A 17.5cm, Aw 16 cm。 育 苗 期 间 ， 
每 株 施用 完全 营养 液 2 次 ， 每 次 10 mL。 当 幼苗 两 叶 一 心 时 ， 选 择 长 
势 ， 大 小 均一 的 幼苗 进行 低温 处 理 。 低 温 胁迫 处 理 为 4“C， 光 照 强 
度 为 5 000 lx， 低 温 处 理 24 h; 对 照 (СК) 温度 为 28“C， 光 照 强 度 
73 5 000 Ix, 
1.2.2 КМА 提取 及 Шшпіпа 测序 

采集 低温 处 理 和 对 照样 本 的 叶片 ， 及 时 用 TRNzol Reagent 试剂 
盒 (天 根 生化 科技 有 限 公 司 , 北京 ) 提 取 并 检测 RNA。 将 低温 胁迫 处 理 
组 与 对 照样 本 组 提取 获得 的 高 质量 RNA 送 往 华 大 基因 科技 有 限 公司 


(深圳 ) 进 行文 库 构 建 ,每 个 处 理 3 个 重复 .采用 Ilumina HiSeq TM 2000 


平台 对 质量 合格 的 RNA 文库 进行 测序 (Ali et al., 2008; Wang et al., 
2009)。 通 过 去 接头 、 去 低 质量 、 去 污染 等 过 程 完 成 数据 处 理 获得 干 
净 序 列 后 ， 使 用 Bowtie2 软件 按照 MiRbase > pirnabank > 
snoRNA (human/plant) > Rfam > other sRNA 的 优先 级 顺序 将 小 RNA 
(SRNA) 遍历 注释 获得 未 注释 RNA HER. 

1.2.3 ВЯ miRNA 的 鉴定 及 新 miRNA 预测 

采用 miRDeep2 软件 对 所 得 的 未 注释 sRNA 序列 与 参考 基因 组 
ГЕНЕ (5ассһағит spontaneum ) 基 因 组 ] 进 行 序列 比 对 分 析 (Marc, 
2008; Maurits et al., 2015), 从 而 鉴定 出 已 知 的 甘蔗 耐寒 相关 miRNAs; 
将 过 滤 获得 的 未 比 对 序列 比 对 至 参考 序列 ， 通 过 碱 基数 目 延 伸 、 
miRNA 结构 预测 的 方法 ， 获 得 新 的 miRNA, 
1.2.4 差异 表达 miRNA 的 鉴定 

对 耐寒 型 甘蔗 样品 的 已 知 miRNA 的 读数 进行 分 析 ， 判 断 低 温 胁 
迫 前 后 不 同 耐寒 型 甘蔗 样品 中 miRNA 的 差异 表达 ， 以 
llog2(FoldChange)>1, P < 0.05 为 筛选 标准 ， 获 得 低温 胁迫 前 后 差异 
表达 miRNA, 
1.2.5 miRNA 靶 基 因 预 测 

采用 psRNATarget(Wu, et al. 2012)、TargetFinder(Fahlgren & 
Carrington 2010) 和 Tapirhybrid(Peer, 2010) 软 件 进行 差异 表达 miRNAs 
靶 基 因 预 测 ， 并 对 预测 靶 基 因 进 GO 数据 库 同 源 性 搜索 ， 确 定 其 参与 
的 信号 传导 及 生物 代谢 途径 。 


1.2.6 荣光 定量 PCR 验证 


选取 14 个 差异 表达 miRNAs RHEA, FJA Premier 5.0 设计 
其 成 熟 miRNA 特异 正 向 引物 .通用 反 向 引物 及 stem loop 5 JS 4), 
然后 利用 LightCycler& 480 Instrument II 进行 RT-PCR 扩 增 ， 设 阴性 
对 照 (不 添加 cDNA 模板 ) 以 监控 可 能 的 污染 。 以 甘蔗 GAPDH 为 内 参 
基因 , 每 个 样品 均 设置 3 次 重复 , 采用 2M 法 计算 基因 相对 表达 。 
2 结果 与 分 析 
2.1 高 通 量 测序 数据 分 析 

选择 3 个 耐寒 性 不 同 甘 蔗 品种 的 幼苗 进行 低温 胁迫 , 每 个 品种 的 
对 照 (CK) 和 低温 处 理 (T) 3 个 重复 , 分 别 取 叶片 共 18 个 样本 进行 高 
通 量 测序 , 构建 低温 胁迫 前 后 sSNA 文库 , 然后 对 原始 测序 数据 进行 
3 端 去 接头 、Trim 低 质 量 、 片 段 大 小 选择 等 质 控 处 理 ， 获 取 高 质量 
Clean Read 序列 ， 结 果 见 表 1。 从 表 1 可 以 看 出 , 在 GR2 的 CK 和 TT 
叶片 中 分 别 挖 扬 出 24 310 558 条 、23 925 673 条 测序 数据 ， 最 终 处 理 
分 别 得 到 21 343 194 条 、21 576 594 条 ， 高 质量 Clean Read 序列 分 别 
占 总 序列 的 87.83%、90.25%， 且 所 构建 的 文库 与 参考 基因 组 进行 比 
对 ， 其 CE 和 T 分 别 有 86.71% 与 84.32% 的 序列 与 参考 基因 组 匹配 ; 
GT28 的 CK 和 T 叶片 中 分 别 挖掘 出 23 915 244 条 、23 528 103 条 测 
序数 据 ， 最 终 处 理 分 别 得 到 22 045 186 条 、21 481 129 条 ， 高 质量 
Clean Read 序列 分 别 占 总 序列 的 92.18% 与 91.28%， 且 所 构建 的 文库 
与 参考 基因 组 进行 比 对 ， 其 CK T 43508 71.75% 与 84.1% 的 序列 
与 参考 基因 组 匹配 ;ROCC22 的 CK 和 TT 叶片 中 分 别 挖 扬 出 27 069 867 


条 、23 631 208 条 测序 数据 , 最 终 处 理 分 别 得 到 21 007 371 条 、20 640 


357 条 ,高 质量 Clean Read 序列 占 总 序列 的 79.5995; 87.35%, 且 所 构 
建 的 文库 与 参考 基因 组 进行 比 对 ,其 CK 和 T 分 别 有 92.62%5 85.67% 
的 序列 与 参考 基因 组 匹配 ( 表 1). 


表 1 样品 测序 数据 统计 


Table 1 Sequencing data statistics of the samples 


E 处 理 样本 编号 总 序列 干净 序列 干净 序列 比 匹配 基因 组 匹配 比例 

Library Treatment Sample ID Total reads Clean reads Percentage (96) Mapped genome Percentage (96) 

GR2 CK GR2 1 24 864 865 21 017 789 84.53 17 335 491 82.48 

CK GR2 2 23 889311 21 068 679 88.19 17 556 964 83.33 

CK GR2 3 24 177 497 21 943 114 90.76 20 697 208 94.32 

均值 Mean 24 310 558 21 343 194 87.83 18 529 887 86.71 

GR2 l 23 612 049 22 167 436 93.88 17 332 696 78.19 

GR2 2 23 331 030 21 905 886 93.89 17 569 575 80.2 

GR2 3 23 294 402 21 590 060 92.68 18 360 373 85.04 

均值 Mean 23 925 673 21 576 594 90.25 18 197 456 84.32 

GT28 CK GI28 1 23 077 345 21 109 698 91.47 15 898 066 75.31 

CK GI28 2 24 699 470 22 556 661 91.32 14 254 489 63.19 

CK GI28 3 23 968 917 22 469 199 93.74 17 247 280 76.76 

均值 Mean 23 915 244 22045186 92.18 15 799 945 71.75 

GI28 1 23 421 689 21 192 257 90.48 16 796 178 79.26 

GI28 2 23 271 839 20 918 101 89.89 18 640 603 89.11 

GI28 3 23 890 781 22 333 031 93.48 18 743 809 83.93 

均值 Mean 23 528 103 21 481 129 91.28 18 060 196 84.1 


Ж 处 理 样本 编号 总 序列 干净 序列 干净 序列 比 匹配 基因 组 匹配 比例 
Library Treatment Sample ID Total reads Clean reads Percentage (%) Mapped genome Percentage (96) 
ROC22 CK ROC22 1 23475518 20335619 86.62 17536453 86.24 
CK ROC22 2 32786885 20072677 61.22 19292598 96.11 
CK ROC22 3 24947198 22613819 90.65 21598010 95.5] 
均值 Mean 27069867 21007371 79.5 19475687 92.62 
ROC22 1 24211511 21189727 87.52 18402946 86.85 
ROC22 2 23491833 20043323 85.32 17201360 85.82 
ROC22 3 23190281 20688022 89.21 17448123 84.34 
均值 Mean 23631208 20640357 87.35 17684143 85.67 


SRNA 长 度 与 不 同 功能 有 关 ， 其 中 21~22 nt 的 主要 与 mRNA Y 
割 和 转录 后 基因 沉默 相关 ，24 nt 的 主要 与 RNA 导向 的 DNA 甲 基 化 
和 转录 基因 沉默 相关 。 我们 对 Total clean reads 进行 统计 ， 发 现 SRNA 
主要 集中 在 21 ~24 nt 之 间 , 但 不 同 抗 塞 能 力 的 材料 之 间 有 一 定 差异 ， 
且 不 同 长 度 的 测序 频率 不 同 (图 1). 
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1 长 度 分 布 图 
Fig. 1 sRNA length distribution 
进一步 对 每 个 重复 样品 之 间 ， 根 据 基因 表达 量 进 行 Pearson 相关 


系数 计算 ， 并 将 这 些 系 数 以 热 图 的 形式 反映 出 来 (图 2)。 从 Pearson 


相关 系数 可 以 看 出 ， 同 一 品种 相同 处 理 的 重复 样品 间 相 关 性 较 高 
(0.667~0.990)， 表 明 所 有 样品 的 表达 量 基本 保持 一 致 ， 说 明 其 符合 
样本 重复 性 实验 标准 ， 可 以 满足 后 续 的 差异 表达 分 析 。 
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X and Y axis represent each sample. The color represents the correlation coefficient. 
The bluer the color, the higher the correlation; and the lighter the color, the lower the 


correlation. 


2 样品 间 相 关 性 分 析 热 图 
Fig. 2 Correlation heatmap ofthe different samples 


2.2 已 知 miRNA 分 析 及 新 miRNA 预测 

利用 AASRA 软件 将 核 苷 酸 序 列 比 对 到 参考 基因 组 及 miRBase 
数据 库 中 , 鉴定 得 到 分 属于 84 个 已 知 miRNA 家 族 的 322 个 miRNA, 
其 中 对 照 组 中 的 297 个 miRNA 分 属于 69 个 家 族 ， 低 温 处 理 组 中 的 
305 个 miRNA 分 属于 74 个 家 族 ,在 这 些 家 族 中 ,拥有 最 多 家 族 成 员 
数量 的 为 miR169(32 个 )， 其 次 分 别 为 miR166、miR171、miR167、 


miR156、miR396 (图 3)。 
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3 miRNA 家 族 数量 分 布 前 18 位 
Fig. 3 The top 18 miRNA family members 
根据 成 熟 植物 miRNA 序列 的 高 度 保守 型 及 miRNA 前 体 拥有 标 
志 性 发 夹 结构 的 特性 ， 利 用 同 源 搜寻 比 对 的 方法 在 割 手 密 miRBase 


数据 库 中 进行 搜索 ， 对 具有 同 源 性 的 miRNA 序列 进行 提取 ， 选 取 其 
中 表达 量 最 高 的 作为 甘蔗 响应 低温 的 保守 miRNA 的 候选 者 ， 最 终 得 
到 甘蔗 响应 低温 的 新 miRNA 110 个 。 
2.3 差异 表达 miRNA 分 析 

过 比较 对 照 组 与 处 理 组 miRNA 的 表达 量变 化 情况 ， 可 以 判断 
低温 胁迫 下 抗 寒 能 力 不 同 的 甘蔗 品种 miRNA 的 差异 表达 情况 。 在 差 
FRA miRNA 检测 过 程 中 , Llllog2(FCO)>2, P<0.05 作为 筛选 标准 ， 
低温 胁迫 时 差异 表达 miRNA 如 表 2 所 示 , 包括 100 个 已 知 miRNAs 
(61 个 上 调 ，39 个 下 调 )，37 个 新 miRNAs (15 个 上 调 ，22 个 下 调 )。 
呈 下 调 趋势 的 miRNA KEH miR8175, miR5564, miR444, miR166 
在 内 的 21 个 家 族 , 呈 上 调 趋势 的 miRNA 家 族 包 括 miR156、miR169、 
miR172、miR393、miR397、miR408 等 10 个 家 族 。 这 些 具有 差异 性 
表达 的 miRNA 在 甘蔗 响应 低温 胁 授 情况 下 可 能 帮 挥 特定 的 功能 。 

表 2 甘蔗 响应 低温 差异 表达 miRNA 


Tab. 2 Differential expression of miRNAs of sugarcane in response to low temperature 


үй 家 族 miRNA 编号 对 照 表达 量 处 理 表 达 量 log2FC P fü 趋 
No. Family miRNA ID Expression (CK) Expression (T) Pvalue Up/Down 
1 miR156 miR156 28. 734 948. 957 5. 045 0.000 35 Up 
miR156c-3p 2 12. 607 112. 513 3. 158 0.000 15 Up 
miR156e-3p 122. 352 13. 266 -3.205 0.000 11 Down 
miRl156f-5p 31. 991 201. 084 2. 652 0.001 10 Up 
miR156k-3p 126. 664 51. 112 -1.309 0.021 58 Down 
2 miR159  miRl59a 1 136 9. 123 620. 576 -1.142 0. 021 27 Down 
miR159a 4 78. 216 15. 693 -2.317 0.000 02 Down 
miR159a-5p 4 241. 182 820 8.194 5. 089 0.000 00 Up 
3 miR160 miR160 332. 983 436 13.722 7. 033 0.000 00 Up 
miR160a-3p 3 68. 014 366. 331 2. 429 0.006 62 Up 
miR160a-3p 4 32. 520 421. 932 3. 698 0.000 00 Up 
miR160e-3p 206. 117 693. 955 1. 751 0.010 93 Up 


x 家 族 miRNA 编号 对 照 表 达 量 处 理 表达 量 log2FC P fü 趋势 
"A" Family miRNA ID Expression (CK) Expression (T) Pvalue Up/Down 
mi R160f-5p 138 0. 969 441. 538 -1.645 0.006 59 Down 
miRl60g 1 185. 534 24. 856 -2.900 0.000 04 Down 
4  miRl62 miR162-5p 231. 769 225 8.160 3.284 0.000 46 Up 
mi R162a-3p 739. 997 28. 331 -4.707 0.001 77 Down 
5 — miRl64 miR164b-3p 2 5. 405 51. 477 3.252 0.000 21 Up 
miR164c-3p 2 4. 473 186. 499 5.382 0.00000 Up 
mi R164£-3p 0. 604 12. 768 4.403 0.00864 Up 
6 — miRl66  miRl166a 365. 383 109. 611 -1.737 0.000 00 Down 
miRl66a-5p 1 35. 859 164. 403 2.197 0.021 08 Up 
miR166d-5p 2 661. 392 566 8.546 3.099 0.000 13 Up 
miRl66h-3p 1 601 9.722 269 4.272 -1.160 0.010 14 Down 
mi RL66k-3p 102 1. 709 469. 918 -1.121 0.004 79 Down 
miRl66m 1 596. 293 124. 197 -2.263 0.000 02 Down 
miR166m 2 204. 395 7. 433 -4.781 0.000 46 Down 
7 — miRl67  miRl67a-5p 654. 093 176. 527 -1.890 0.000 20 Down 
miR167d 540 2. 076 579 07.162 3.422 0.000 14 Up 
miRl67d 1 190. 952 21.143 -5.207 0.000 00 Down 
miRl67d-3p 5 4. 235 37. 787 3.157 0.003 92 Up 
mi R167g-3p 4. 047 21. 062 2.380 0.003 16 Up 
8 — miRi68 miRl68a-3p 2 1749 8.058 362 33.268 2.273 0.000 03 Up 
9 шіЮІ60  miRl69a-3p 3 1. 682 36. 573 4.443 0.000 74 Up 
miRl69c-3p 1 10. 808 65. 672 2.603 0.000 17 Up 
miRl69c-3p 5 25.511 579. 499 4.506 0.000 00 Up 
miR169d-3p 3 357. 242 461 8.614 3.692 0.000 28 Up 
miRl69d-5p 1 65. 463 107 9.187 4.043 0.001 27 Up 
miRl69e 2 20. 181 2. 225 -3.181 0.003 14 Down 
miR169e 3 3. 082 26. 689 3.114 0.004 26 Up 
miRl69e-3p 1 58. 233 241. 058 2.049 0.010 57 Up 
miR169h 2 7.341 93. 031 3.664 0.006 36 Up 
niRl69k-5p 161. 245 721. 733 2.174 0.001 56 Up 
mi R169m 3. 793 76. 157 4.328 0.000 14 Up 
mi Rl69m-5p 1.199 18. 117 3.917 0.002 27 Up 
mi R169r-3p 3. 887 78. 355 4.333 0.000 04 Up 
miR169r-5p 2. 751 68. 185 4.631 0.000 00 Up 
10 шіРІТІ miRl7la 3 65. 519 23. 729 -1.465 0.017 88 Down 
miR171a-3p 15. 331 4. 377 -1.808 0.020 04 Down 
miR171b-3p 3 137. 190 170 3. 706 3.628 0.003 49 Up 
miRl7le-5p 1 9. 596 43. 006 2.164 0.006 55 Up 
miR171f 8.261 0.281 -4.879 0.001 25 Down 
miR171i 1 222 8. 766 818. 757 -1.445 0.000 02 Down 


Т 家 族 miRNA 编号 对 照 表 达 量 处 理 表达 量 P fü 趋势 
Family miRNA ID Expression (CK) Expression (T) s Pvalue Up/Down 
11 miR172 miRl72a 2 12. 532 1. 849 -2.761 0.010 41 Down 
miR172b-5p З 346. 643 276 3. 508 2. 995 0.001 70 Up 
miR172d-5p_4 105. 771 523. 500 2. 307 0.013 72 Up 
12  miR2118 miR2118a 2 26. 740 4. 574 -2.548 0.001 31 Down 
miR2118b 533 9. 363 209 7.461 -1.348 0.014 70 Down 
miR2118d 96. 120 402. 856 2. 067 0.000 81 Up 
miR2118f 1 90. 692 521. 181 2. 523 0.006 67 Up 
miR2118g 66. 283 5. 469 -3.599 0.000 01 Down 
miR2118p 498. 943 1827. 679 1. 873 0.012 49 Up 
13  miR390 miR390-3p 2 2. 069 15. 562 2.911 0.004 43 Up 
14 miR393 . miR393-3p 1 53. 725 698. 525 3. 701 0.000 01 Up 
miR393a 3 107. 640 542. 872 2. 334 0.000 34 Up 
15  miR395 miR395h-5p 2. 586 30. 738 3. 571 0.008 14 Up 
16 miR396 шіК396а-3р 4 1. 601 25. 432 3. 989 0.003 31 Up 
miR396a-5p 180. 860 17. 171 -3.397 0.000 51 Down 
miR396b 13. 610 271. 022 4. 316 0.000 00 Up 
miR396b-3p 1 6. 136 13. 003 3. 573 0.000 46 Up 
miR396b-3p 3 35. 226 171. 642 2. 334 0.000 78 Up 
miR396b-5p 183. 564 51. 783 -1.826 0.000 80 Down 
miR396e-3p 1 182. 946 44. 050 -2.054 0.010 03 Down 
miR396g 138. 206 55. 367 -1.320 0.010 45 Down 
miR396h 7.617 1. 352 -2.495 0.012 08 Down 
17 miR397 . miR397-3p 3 215. 647 263 7. 398 3.612 0.001 77 Up 
miR397a 3 27. 767 336. 211 3. 598 0.001 87 Up 
18 miR398 | miR398a-3p 2 11. 259 119. 098 3. 403 0.005 51 Up 
miR398b 20. 516 160 8.552 6. 293 0.000 03 Up 
19 miR399  miR399 1 16. 048 1. 571 -3.353 0.003 10 Down 
miR399b-5p 1 14. 959 76. 405 2. 353 0.010 00 Up 
miR399e-5p 19. 589 295. 499 1. 893 0.002 30 Up 
miR399k 1 35. 366 154. 991 2. 132 0.000 73 Up 
20 miR408 піК408-3р 1 951. 933 111 95. 030 3. 556 0.020 29 Up 
miR408-5p 9 6. 759 599. 555 6. 471 0.000 00 Up 
miR408b 1 16. 628 467. 667 4. 814 0.000 39 Up 
mi RA08d 116. 205 230 6.584 4. 311 0.000 71 Up 
2] | miR444 miR444b. 2 946 6. 966 226 9.188 -2.061 0. 001 85 Down 
22 miR4995  miR4995 297. 645 71. 203 -2.064 0.018 57 Down 
23  miR5139 miR5139 43. 723 7. 385 -2.566 0.010 29 Down 
24 miR529 miR529-3p 42. 770 5. 191 -3.042 0. 000 55 Down 
25  miR530 miR530a 2 24. 562 4. 006 -2.616 0.017 36 Down 
26 | miR5505  miR5505 13. 384 0. 482 -4.795 0.000 04 Down 
27  miR5523 | miR5523 9.527 0.604 -3.979 0.002 70 Down 
28 miR5564 miR5564c-5p 205 7.141 274. 385 -2.906 0.000 00 Down 
29  miR5568  miR5568c-3p 2. 390 23. 041 3. 269 0.003 51 Up 
30  miR6218 | miR6218-5p 4.570 20.811 2.187 0.008 34 Up 
31  miR6221 miR6221-5p 50. 052 12. 539 -1.997 0.009 55 Down 
32  miR6223  |miR6223-5p 14. 898 249. 912 4. 068 0.000 00 Up 
33  miR6224 miR6234a-3p 0. 548 7. 645 3. 801 0.004 88 Up 
34  miR8175 miR8175 707. 142 270. 620 -1.386 0.007 87 Down 


- 家 族 miRNA 编号 对 照 表 达 量 处 理 表达 量 log2FC P fH 趋势 
Ns Family miRNA ID Expression (CK) Expression (T) Pvalue Up/Down 
35 n miri04 novel mir104 616 0.659 415. 650 -3.890 0.000 05 Down 
36 n mirl06 novel mir106 4T. A11 1153. 668 4. 605 0.000 02 Up 
37 n mirl07 novel mirl07 9. 744 171. 505 4. 138 0.000 24 Up 
38 n miri08 novel mir108 891. 184 72. 596 -3.618 0. 000 09 Down 
39 n пігі09 novel mir109 16. 716 194. 103 3. 538 0.001 33 Up 
40 n miril novel mirll 660. 280 325 8.241 2. 303 0.003 16 Up 
4] n mirll0 novel тіг110 21.128 0. 344 -6.302 0.000 01 Down 
42 n mirl5 novel mirl5 17.837 0. 338 -5.724 0. 000 08 Down 
43 n mirl7 novel mirl7 11. 582 0. 296 -5.288 0.000 33 Down 
44 n mir20 novel mir20 0. 748 142. 501 1. 573 0.000 00 Up 
45 n mir23 novel mir23 23. 022 0. 352 -6.032 0.000 02 Down 
46 n mir24 novel mir24 248. 551 14. 516 -1.738 0.012 12 Down 
47 n mir25 novel mir25 241. 891 188. 583 1. 705 0.003 72 Up 
48 n mir26 novel mir26 148. 771 43. 014 -1.790 0. 001 28 Down 
49 n mir27 novel mir27 11. 000 0. 643 -4.096 0.008 96 Down 
50 n mir29 novel mir29 204 4. 588 335. 897 -2.606 0.000 34 Down 
51 n mir3 novel mir3 12. 479 0. 790 -3.982 0.001 83 Down 
52 n mir34 novel mir34 0. 545 20. 981 5.267 0.000 89 Up 
53 n mir36 novel mir36 37. 565 0. 505 -6.218 0.000 01 Down 
54 n mir37 novel mir37 20. 896 662. 275 4. 986 0.00000 Up 
55 n mir42 novel mir42 102 2.840 358. 059 -1.514 0.005 06 Down 
56 n mir47 novel mir47 13. 512 0. 680 -4.313 0.013 92 Down 
57 n mirbl novel mirb5l 0. 535 121. 564 1.828 0.000 00 Up 
58 n mir53 novel mir53 577. 260 169. 427 -1.769 0. 009 58 Down 
59 n mirb8 novel mir58 68. 047 926. 881 3. 768 0.006 66 Up 
60 n mir61 novel mir61 191 2.408 432. 242 -2.145 0.000 00 Down 
61 n mir62 novel mir62 180. 637 63. 979 -1.497 0. 020 85 Down 
62 n mir64 novel mir64 7. 505 0. 474 -3.984 0. 024 09 Down 
63 n mir65 novel mir65 0. 618 37. 647 5. 929 0.000 00 Up 
64 n mir72 novel mir72 0. 828 260. 425 8. 296 0.000 00 Up 
65 n mir75 novel mir75 186. 658 37. 721 -2.307 0.001 83 Down 
66 n mir82 novel mir82 190. 143 854. 423 2. 168 0.011 08 Up 
67 n mir84 novel mir84 12. 396 0. 662 -4.221 0. 006 72 Down 
68 n mir86 novel mir86 0. 439 11. 476 4. 709 0.003 54 Up 
69 n mir89 novel mir89 5. 958 0. 419 -8.828 0.021 37 Down 
70 n mir95 novel mir95 294. 205 33. 471 -3.136 0.008 67 Down 
71 n mir96 novel mir96 0. 862 51. 926 5. 912 0.000 16 Up 
2.4 miRNA (ЕВЕ 90155917 

根据 筛选 差异 表达 的 miRNA 与 对 应 物种 的 基因 序列 信息 ， 使 用 


psRobot, TargetFinder 和 Tapurhybrid RPH TERA, 结果 见 图 
4, 3 种 预测 方法 共 预 测 了 1 844 个 相同 的 潜在 靶 基 因 ， 其 中 1 696 个 
属于 已 知 miRNA 靶 基 因 ，148 个 属于 新 miRNA 靶 基 因 。 为 进一步 分 
析 预 测 所 得 靶 基 因 的 潜在 生物 学 功能 ， 对 预测 所 得 的 靶 基 因 进行 GO 
分 析 ， 共 鉴定 出 甘蔗 响应 低温 胁迫 miRNA 靶 基 因 在 生物 过 程 

(Biological process) 有 13 个 功能 亚 类 ， 在 细胞 组 分 (Cellular 
component) 有 11 个 功能 亚 类 ， 在 分 子 功 能 (Molecular function) 有 
7 个 功能 亚 类 (图 4)。 在 生物 过 程 类 别 中， 预测 靶 基因 的 主要 生物 学 
功能 富 集 于 细胞 过 程 和 代谢 过 程 ; 在 细胞 组 分 类 别 中 , 预测 靶 基 因 生 
物 学 功能 主要 语 集 于 细胞 、 细 胞 器 和 生物 膜 ; 在 分 子 功 能 类 别 中 ,， 预 
测 计 基因 生物 学 功能 主要 语 集 于 结合 和 众 化 活性 及 转运 活性 。 绝 大 多 
数 的 靶 基 因 功 能 均 与 这 些 结合 功能 及 其 他 相近 的 结合 功能 相关 。 大 多 
数 miRNA 通过 直接 或 者 间接 介 导 靶 基 因 的 表达 调控 相关 代谢 途径 响 
应 低温 胁迫 ， 这 些 miRNA 所 调控 的 靶 基 因 对 甘蔗 的 耐寒 性 起 着 关键 
的 调控 作用 ( 表 3)。 


mA HE m 细胞 组 分 分 子 功能 
Biological process Cellular component Molecular function 


基因 数量 Number of genes 


1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 


4 差异 表达 miRNAs 靶 基 因 的 GO 功能 分 析 
Fig. 4 GO analysis of target genes for differentially expressed miRNAs 
Т. 细胞 过 程 ; 2. 代谢 过 程 ;3. 生 物 调 节 ; 4. 发 育 过 程 ; 5. 定位 ; 6. 响应 刺激 ; 7. 多 


细胞 有 机 体 过 程 ; 8. 繁殖 ; 9. 繁殖 过 程 ; 10. 信号 传导 ; 11. 细胞 成 分 组 织 或 生物 

RE; 12. 多 有 机 体 过 程 ; 13. 生长 ; 14. 细胞 Cell; 15. 细胞 器 ; 16. 细胞 膜 ; 17. ZH 

胞 膜 部 分 ; 18. 蛋白 复合 体 ; 19. 细胞 器 部 分 ; 20. 细胞 连接 ; 21. 共 质 体 ; 22. 膜 包 

管 腔 ; 23. 超 分 子 复合 体 ; 24. 胞 外 区 ; 25. 结合 ; 26. 转录 因子 活性 ; 27. 催化 活性 ; 
28. 转运 活性 ; 29. 结构 分 子 活性 ; 30. 分 子 功 能 因子 ; 31. 抗 氧 化 活性 . 

1. Cellular process; 2. Metabolic process; 3. Biological regulation; 4. Developmental 


process; 5. Localization; 6. Response to stimulus; 7. Multicellular organismal process; 


8. Reproduction; 9. Reproductive process; 10. Signaling; 11. Cellular component 
organization or biogenesis; 12. Multi-organism process; 13. Growth; 14. Cell; 15. 
Organelle; 16. Membrane; 17. Membrane part; 18. Protein-containing complex; 19. 
Organelle part; 20. Cell Junction; 21. Symplast; 22. Membrane-enclosed lumen; 23. 
Supramolecular complex; 24. Extracellular region; 25. Binding; 26. Transcription 
regulator activity; 27. Catalytic activity; 28. Transporter activity; 29. Structural 


molecule activity; 30. Molecular function regulator; 31. Antioxidant activity. 


353 部 分 甘蔗 响应 低温 miRNAs 及 其 靶 基 因 


Tab. 3 Some potential target genes for miRNAs in sugarcane responding to low temperature 


miRNA 编号 “表达 趋势 Зо ТАС; 期 望 值 功能 描述 
miRNA ID Expression Targets ID Expectation Functional description 
T fle» aH 启 云 结合 13 
miR156 Nus Sspon.03G0000360-3C 3 оа 


Squamosa promoter-binding-like protein 13 


生长 素 啊 应 因子 22 
Sspon.01G0019580-1A 2 ñ 长 素 有 应 因子 
Auxin response factor 22 
! 下 调 生长 素 响应 因子 8 
miR160g_1 Down Sspon.04G0008020-4D 3 - Кж Шен 
Auxin response factor 8 
前 11212 ` HI 1 , П % 
Sspon.01G0032360-1A 3 Ы жын тА . 
Prolycopene isomerase 1, chloroplastic 
Wi 263) 3 Hg lm 型 ХІ 
miR167d 2. Sspon.01G0032360-2D 3 RRE PRA " Fi 
p Carotenoid isomerase 1 isoform X1 
KAT Y ЖЖ А-7 ХІ 
Sspon.01G0004210-1A 2:5 FRET a 异 构 体 
Nuclear transcription factor Y subunit A-7 isoform X1 
调 KAF Y T А-4 
miR169a-3p_3 s Sspon.01G0040290-2D 2.5 d А 5 
р Nuclear transcription factor Y subunit A-4 
录 因 子 Y 亚 基 A-3 
Sspon.02G0034130-2C 2.5 apad жо 
Nuclear transcription factor Y subunit А-3 
ü ЦЕ 
、 Sspon.01G0000260-3P 3 ЛЕ н : 
miR171b-3p 3 上 调 Scarecrow-like protein 6 
Б 26 Ор EE ИЖ 27 
Sspon.04G0006810-1A 3 РЕЛШЕ 


Scarecrow-like protein 27 
miR172d-5p 4 上 调 Sspon.02G0011120-1A 3.5 转录 因子 bHLH49 


miRNA 编号 “表达 趋势 RUE DAR s 期 望 值 功能 描述 
miRNA ID Expression Targets ID Expectation Functional description 
Up Transcription factor bHLH49-like 
RNA fi 40A 
$spon.03G0021730-1A 3.5 Еа Шеик. 
Pre-mRNA--processing protein 40A 
Ë /结合 蛋 E 
Sspon.04G0008430-4D 3.5 а СА 
microtubule-associated protein futsch 
DUF1682 ЖИЕН 
S .07G0023470-1B 3 
SESS DUF1682 family protein 
ys Ж П ЧИ p-4 
miR2118d 2 Sspon.01G0014780-1A 3 НЕНИЯ Ps : 
P Putative casein kinase II subunit beta-4 
: 上 调 а -L- 阿 拉 伯 糖苷 酶 1 
R393-3p_1 Е 20400-ІР 
Е Ор Se i x Alpha-L-arabinofuranosidase 1 
ji 合 系统 ПАЖ E 2 
miR396b 4. Sspon.04G0022140-2C 3 а . ; 
P Photosystem II oxygen-evolving enhancer protein 2 
RIE SAG39 
f Sspon.01G0027900-2C 3 EERE B T ЕГА 
;R397a 3 上 调 Senescence-specific cysteine protease SAG39 
mi a 
= Up 浓缩 素 -2 复合 物 亚 基 рз 
Sspon.04G0017800-3D 3 
Condensin-2 complex subunit D3 
上 调 硒 结合 НЕҢ ЕН 1 
IR. : 28270-2P 3 
шы Up — MIS Selenium-binding protein 1 
上 调 泛 素 结合 酶 E2 23 
iR399k 1 ; 21800-ІВ 3.5 
а Ор ааа Probable ubiquitin-conjugating enzyme Е2 23 
类 质 体 蓝 素 前 体 
Sspon.01G0047030-3D 3.5 joi 
;«RA08d 上 调 Chemocyanin precursor 
mi 
Up 小 亚 基 核糖 体 S23 
Sspon.07G0017290-2B 3 


Small subunit ribosomal protein S23 


miRNA 编号 ”表达 趋势 RUE DAR s 期 望 值 功能 描述 
miRNA ID Expression Targets ID Expectation Functional description 
下 调 Sspon.03G0002860-2B 3.5 核糖 核酸 酶 I， 叶 绿 体 /线粒体 异 构 体 ХІ 
miR5140 D ies Ribonuclease II, chloroplastic/mitochondrial isoform X1 
m 
Sspon.07G0015510-1A 3 add f f ! 
Transposon protein, putative, unclassified 
凋 亡 前 期 丝氨酸 和 蛋白酶 NMA111 
f Sspon.02G0012470-1A ^ 3.5 ОРЕВ 
Raira 下 调 Pro-apoptotic serine protease ММА 111 
mi a Lib 
I “WD Sspon.05G0000530-2D 35 CERA 结合 域 的 蛋白 质 1 
pons ` Acyl-CoA-binding domain-containing protein 1 
l 下 调 泛 素 羧基 端 水 解 酶 
R5523 .02 1670-3D 3.5 
Down са ды Ubiquitin carboxyl-terminal hydrolase 
. 下 调 3'-5- 核 糖 核酸 外 切 酶 家 组 蛋 F 
R5564c-5 Р 2720-30 3 
= Р Down аза 3'-5'-exoribonuclease family protein 
ЖЕЛДЕЙ ХІ 
Sspon.01G0013100-1A 3 HE 
Probable trehalase isoform X1 
二 酰 辅酶 A 酰基 载体 蛋 
Sspon.01G0014700-3C 3 PM ВЕ ЕЕЕ 
Malonyl CoA-acyl carrier protein transacylase 
5- FER CR 
f 上 调 Sspon.01G0033120-1A 2.5 ыы 
miR5568c-3p Up 5-oxoprolinase 
贡 网 相关 的 降解 蛋 
Sspon.01G0033150-1P 3.5 шана аа 
Derlin-2.2 
Sspon.03G0010210-1P 3 RRHANOREUNEUM 2 2 f 
Malate dehydrogenase 2 mitochondrial 
Sspon.04G0006870-1A 3.5 pe WB ERAI 1 


miRNA 编号 “表达 趋势 RUE DAR s 期 望 值 功能 描述 
miRNA ID Expression Targets ID Expectation Functional description 
Probable trehalose-phosphate phosphatase 1 
ABC 转运 蛋 F 族 成 员 22 同型 X2 
、 Sspon.01G0003280-1A 3.5 iini iid n : 
Е 上 调 ABC transporter G family member 22 isoform X2 
s Up 605 核糖 体 蛋 白 136 
Sspon.02G0047980-1C 3.5 diis A . : . 
605 ribosomal protein L36a, partial 
i 过 氧 酶 体 生 成 蛋白 -14 
miRG234a3p $ ^ Sspon07G0037700-ID з ЛН 
р Peroxin-14 
AE H M TE 
Sspon.01G0023790-2P 3.5 EKR AWRA А 
Auxin-responsive protein IAA30 
DNA 解 旋 酶 IN 
Sspon.01G0039360-1P 3.5 D dí pi 
DNA helicase ІМО80 
ПЕ PUEA КМ51 
novel mir23 Ki B Sspon.02G0004350-3D 3.5 НОШКЕН . 
own Vacuole membrane protein KMS1 
ЕЗ 泛 素 蛋白 连 AIP2 
Sspon.02G0014140-3D 2 E ue 
E3 ubiquitin-protein ligase AIP2 
不 乳 碱 性 亮 氨 酸 拉链 转录 激活 因子 5 
Sspon.02G0014480-3C 5 ЛЕЕ P жен ALT 
Opaque endosperm 5 
К» х 33. pa = zL E 
l Sspon.04G0017230-1A 2.5 SHROSHSHTÉRGEDS | f 
—À TT US 4Squamosa promoter-binding-like protein 4 
- Down 50S 核糖 体 和 蛋白 133 
Sspon.07G0022740-1B 3.5 RE 
508 ribosomal protein L33 
iji Ж O- 和 S- 糖 基 复 合 
novel mir90 Lu Sspon.08G0003790-2B 2.5 Е О- Ж SCHGER EI 
- Down Hydrolyse O- and S-glycosyl compounds 


2.6 qRT-PCR 验证 

利用 qRT-PCR 技术 对 测序 所 得 到 的 miRNA 及 靶 基 因 的 差异 丰 度 
进行 验证 ， 而 miRNA 及 其 靶 基 因 的 引物 序列 如 表 4 所 示 。 在 对 照 组 
与 低温 处 理 中 共 筛 选 出 14 个 差异 表达 的 miRNA 及 其 靶 基 因 对 其 进 
{T RI-PCR 验证 ,包括 miR156.miR160g 1, miR167d, miR169 a-3p 3, 
miR171b-3p 3, miR172d-5p 4, miR393-3p 1, miR396b, miR397a 3, 
miR398b, miR399k 1, miR408d, novel mir36, novel тіг89( 5), 
根据 表 4 和 图 5 的 结果 , 除 novel-miR36 外 ,选择 的 其 余 13 个 miRNA 
在 qRT-PCR 实验 中 的 表达 模式 与 高 通 量 测序 检测 到 的 表达 模式 一 致 
(图 5: A), 这 表明 大 部 分 高 通 量 测序 的 结果 均 能 被 艾 光 定量 PCR ix 
术 所 验证 , 此 次 研究 数据 为 真实 可 信和 的。 所 筛选 的 14 个 miRNA rh, 
BR novel-miR36 ДЬ, 其余 13 个 miRNA 均 与 其 靶 基 因 呈 负 调 控 关 系 
(图 5: В), ШЕҢ miRNA 通常 以 负 调 控 或 沉默 靶 基 因 的 方式 来 调控 植 
物 的 生长 、 发 育 及 对 环境 的 应 激 反 应 。 


424 miRNAs 和 靶 基 因 qRT-PCR 引物 设计 
Tab. 4 The primer of miRNAs and target's gene for qRT-PCR 


引物 
miRNA 编 引物 靶 基 因 编 号 primer 
号 ИИ Target genes 正 向 引物 (5 -3 ) 反 向 引物 (5 -3”) 
miRNA ID ID Forward primer Reverse primer 
(5 3 а= 3 
aise CGCTGACAGAAGAGAGT Sspon. 03G0000 GGCACAACTGAGAGC ACTACATTGGATGGC 
GAGCAC 360-3C ACACCTT AGCACCT 
miRl60g . TATATGCCTGGCTCCTT 5ѕроп. 0160019 | GGAGCACCTTCGTCA АТССАСТССАТСССА 
1 GTATGCCA 580-1A ACCACAA CCACAG 
: TGAAGCTGCCAGCATGA — Sspon. 0160032 ССАССАТССТТССТТ | TGACACCAACTGCTC 
ае TCTGG 360-2D GAGTACC GACCATT 
miR169a- ССАТАСССААССАТСАС — Sspon. 0160004 — ACGAGCGAAGCAAAG АСАССААСТСАСССА 
3p 3 TTGCCG 210-1A CAAGGT GGAAGGA 
miRl71lb- | CTTGAGCCGTGCCAATA | Sspon.01G0000 — TGGATCATTGGCGGC САССААСТТСССССТ 
3p 3 TCACG 260-3P GAGGA GTGGA 
miR172d- ССССАССАССАТСААСА 5<роп. 0260011 | GATGGTGCTGGTGTC ССТССТСТАССТТСС 
bp 4 TTCAC 120-1А GTTCTGA GCTTCAC 
miR393-3 ССССАСТССААТСССТТ 5ѕроп. 0560020 | GCTTCGCCTCCTGCT ТССАССАТСССАТАС 
р 1 ТССААТТ 400-1P TCTTCT CTCCTCT 
. CGCGTTCCACAGCTTTC — Sspon. 0460022 — AGGGTTACACGATTC | CGCCTTCTGCTCTGC 
DP TTGAACT 140-2C CGCTTGA CACAT 
miR397a |^ CGTCATTGAGTGCAGCG 5ѕроп. 0460017 | COGCCAATAAGGTCG — GGTACGCCTTGTGCA 
3 TTGATG 800-3D TCGTCAC GCTTCTC 
niR398b TATATAGGGGCGGACTG 5ѕроп. 0360028 | CGCCGCTTCCTGATC ССАСАТАСАСАССАС 
GGAACAC 270=2P TTACCTT CTGACCT 
miR399k TGCCAAAGGAAATTTGC Sspon.07G0021 TGGAGGAATTGGAGC АСАССССАССАСССА 
1 CCCG 800-1B CGTTGGA AGAAGAG 
niR408d TGCACTGCCTCTTCCCT — Sspon. 0760017 — GATGCCGTTAGGGTC АССАСАССССТТССА 
GG 290-2B AGGACAA ACACAA 
novel mi |GCGATTGACAGAAGAGA Sspon. 0860015  TCGCTTGGATCGTTG АСАТАСССТССТССТ 
r36 GTGAGCAC 870-2В AGAATTG AGTACC 
novel mi TCGACTCCCACTGTGGT 5<роп. 0860003 TCGGCGTCCCTGGTG САСАСССТСААССТС 
r89 CG 790-2В ТТСАА TCGTTCC 


cT {ii 


CT value 
CT 全 
CT value 


A. miRNA qRT-PCR 检验 ; B. 靶 基 因 qRT-PCR 检验 。 
A. The results of miRNA qRT-PCR; B. The results of target genes. 


5 miRNA 及 其 靶 基 因 qRT-PCR 分 析 
Fig. 5 qRT-PCR analysis ofmiRNA and target genes 


3 讨论 与 结论 
miRNA 是 在 植物 中 广泛 分 布 的 一 类 非 编码 小 RNA 分 子 , 在 植物 
生长 、 友 育 和 对 环境 应 激 反应 中 上 友 挥 着 关键 的 调控 作用 ( 刘 运 华 等 ， 


培 品种 以 热带 种 的 种 质 为 主体 。 作 为 广西 传统 的 主导 产业 , 而 低温 不 
仅 是 限制 其 扩大 种 植 区 域 和 实现 高 产 稳 产 的 重要 因素 之 一 , 还 影响 蔗 
农 收益 和 糖 业 稳定 发 展 ( 苏 永 秀 等 , 20060)。 为 了 解 甘蔗 响应 低温 的 内 
在 分 子 机 制 ， 挖 掘 其 与 耐寒 相关 的 miRNAs 及 相关 靶 基 因 ,， 本 研究 对 
不 同 基因 型 甘蔗 进行 低温 胁迫 处 理 , 通过 高 通 量 测序 技术 及 生物 信息 
学 方法 ， 系 统 分 析 不 同 抗 寒 能 力 的 甘蔗 对 低温 胁 据 的 响 。 

通过 高 通 量 测序 技术 及 生物 信息 学 分 析 , 我 们 发 现 SRNA 主要 集 


中 在 21-24 nt Z |8], 但 不 同 抗 寒 能 力 的 材料 之 间 有 一 定 差 寞 ， 且 不 同 
长 度 的 的 测序 频率 不 同 。 许 多 研究 表明 ，sRNA 在 不 同 物种 间 的 长 度 


分 布 会 有 所 区 别 ， 如 拟 南 芥 (Pasquinelli et al., 2000)、 小 麦 (Meng et al., 
2013), 棉花 (Sripathi et al., 2014) 的 SRNA 长 度 最 多 分 布 在 24 nt, 杨 树 
( 李 明 娜 等 , 2014)、 大 豆 (Turner et al., 2012)FTIESAR(Pilcher et al., 2007) 
的 SRNA 长 度 最 多 分 布 在 21 nt， 这 些 结果 与 我 们 的 研究 结果 高 度 相 

似 。 我 们 通过 生物 信息 学 手段 共 挖掘 137 个 miRNA 在 低温 胁迫 前 后 
进行 差异 性 表达 ， 其 中 100 个 已 知 miRNA (61 个 上 调 ，39 个 下 调 )， 
37 个 新 miRNA(15 个 上 调 , 22 个 下 调 )。 许 多 研究 表明 ， 在 低温 胁迫 

下 植物 通过 调节 miRNA 的 表达 水 平 ， 进 而 调节 对 应 靶 基 因 的 表达 ， 
从 而 引起 相关 代谢 与 信号 转 导 途 径 的 变化 来 实现 对 逆境 的 响应 , 其 信 
号 转 导 途径 主要 包括 胞 外 信号 途径 、 胞 内 第 二 信使 、 转 录 因 子 以 及 功 
能 基因 等 。Wu & Poethig (2006) 的 研究 友 现 ， 低 温 胁迫 导致 miR156 

下 调 而 其 靶 基 因 表达 量 增 加 , 进而 调控 拟 南 芥 曹 养 期 延长 、 生 长 代谢 
变 缓 ， 以 应 对 不 良 环境 ， 这 与 本 研究 低温 胁迫 后 miR156 下 调 的 结果 
一 致 。 王 丽 丽 等 (2017) 分 析 在 拟 南 芥 、 小 麦 和 水 稻 研 究 中 发 现 miR160 
响应 低温 胁迫 。 本 研究 中 ,在 低温 胁迫 后 miR160g-1 RATA, WE 
因 表 达 受 抑制 且 作 用 于 生长 素 信号 通路 从 而 在 抵抗 低温 中 发 挥 作 用 。 
Pourcel et al. ( 2005) 研 究 表明 miR397 НЫНАН ЕН, 与 细胞 壁 
木质 素 的 合成 、 抗 病 、 对 环境 的 适应 过 程 密切 相关 。 而 我 们 的 研究 发 
现 低温 胁迫 后 甘蔗 叶片 miR397a-3 呈 上 调 表达 , 其 对 基因 与 抗坏血酸 
氧化 酶 (L-ascorbate oxidase) 有 关 ， 可 能 参与 调控 抗坏血酸 氧化 酶 的 


表达 而 增强 低温 响应 能 力 。 Li et al. 014) 255r miRNAs 对 低温 胁 
迫 的 响应 ,发 现 miR160, miR164, miR394, miR395 和 miR408 具有 
表达 差异 。Sunkar & Zhu (2004) 在 研究 分 析 拟 南 芥 响应 胁迫 miRNAs 
H, 发 现 一 些 miRNA 能 被 多 种 胁迫 因素 诱导 , 如 miR393 受到 低温 、 
干旱 、 高 盐 的 诱导 表达 。Gupta et al. (2014) 研 究 小 麦 mRNAs 对 低温 
PNE. ARE, BANSKI, 发 现在 盐 肋 授 和 低温 胁迫 下 , miR168、 
miR397 均 表 达 下 调 ， 而 miR172 表达 上 调 ; miR393 在 渗透 和 直 胁 迫 
下 表达 量 上 升 ， 在 低温 胁迫 下 表达 量 下 降 。Sun et al.( 2015) 研究 发 
现 ， 在 低温 胁迫 下 对 葡萄 的 miR169 的 表达 量 上 调 ， 但 在 拟 南 芥 、 扁 
桃 中 研究 发 现 miR169 的 表达 量 下 调 。 结 合 前 人 的 研究 结果 ， 我 们 认 
为 特定 miRNAs 对 低温 胁迫 的 响应 可 能 因 植 物种 类 、 同一 植物 的 不 同 
基因 型 、 不 同 的 组 织 类 型 、 胁 迫 的 时 间 等 而 有 所 差异 。 

本 研究 利用 psRNATarget, TargetFinder 及 Tapirhybrid 三 种 分 析 
方式 预测 靶 基 因 ， 对 预测 所 得 的 靶 基 因 进 行 GO 分 析 ， 发 现在 生物 过 
程 类 别 中 预测 靶 基 因 的 主要 生物 学 功能 富 集 于 细胞 过 程 和 代谢 过 程 ; 
在 细胞 组 分 类 别 中 预测 靶 基 因 生 物 学 功能 主要 富 集 于 细胞 、 细胞 器 和 
生物 膜 ; 在 分 子 功 能 类 别 中 预测 对 基因 生物 学 功能 主要 富 集 于 结合 和 
催化 活性 及 转运 活性 。 在 逆境 胁迫 条 件 下 植物 通过 调节 miRNA 的 表 
AKE, 进而 调节 对 应 靶 基 因 的 表达 ,从 而 引起 相关 代谢 与 信号 转 导 
途径 的 变化 来 实现 对 逆境 的 响应 。 因 此 , 进一步 对 参与 调控 植物 激素 
言 号 传导 、 光 合 色 素 合成 、 搞 氧化 酶 系统 、 泛 素 介 导 蛋 晶 水解 、 演 粉 
与 蔗糖 代谢 等 通路 的 靶 基 因 进 行 qRT-PCR 检验 。 本 研究 发 现 , miR156 


靶 向 调控 SBP 转录 因子 ， 低 温 胁迫 前 后 ，miR156 下 调 表达 ， 负 向 调 
控 SBP 转录 因子 ,使 甘蔗 生长 代谢 变 缓 而 耐寒 能 力 增 强 。 梅 琳 (2007) 
研究 冬小麦 发 现 ， 低 温 胁迫 下 miR160 的 靶 基因 ARF17 在 过 表达 
miR160f 拟 南 芥 植株 中 表达 量 也 明显 下 降 。 在 拟 南 芥 、 水 稻 、 玉 米 研 
究 中 也 证 明 miR160.miR167 的 靶 基 因 是 4RF(Auxin response factors) , 
主要 通过 调控 生长 素 信号 通路 应 对 逆境 。 而 我 们 的 研究 也 发 现 甘蔗 在 
低温 胁迫 后 miR160g-1 上 调 表达 且 其 作用 于 转录 因子 ARF。 张 译 云 
等 2012) 研 究 发 现 毛 白杨 在 低温 胁迫 下 miR169ac 表达 下 调 ， 对 其 靶 
基因 转录 因子 NAC 负 向 调控 。 党 春 艳 (2013) 在 研究 局 山 离子 芥 低温 
胁迫 响应 中 ， 发现 miR169a 在 冷 胁迫 下 无 显著 表达 ， 而 miR169 在 冷 
胁迫 下 表达 下 调 ， 表 明 miR169 家 族 之 间 响 应 低温 存在 差异 。 本 研究 
中 miR169 a-3p 3 在 低温 胁迫 后 表达 上 调 ， 作 用 于 细胞 核 转录 因子 使 
其 表达 下 调 ,推测 miR169 а-3р 3 与 调控 甘蔗 低温 耐 受 能 力 密切 相关 。 
此 外 ， 本 研究 通过 生物 信息 学 挖掘 甘蔗 响应 低温 miRNAs， 不 仅 
参与 抗 氧 化 酶 系统 、 植 物 激素 信号 传导 、 遗 传 信息 等 ， 还 发 现 参与 类 
胡 枕 卜 素 代谢 、 中 啉 与 叶绿素 代谢 、 泥 粉 与 蔗糖 代谢 途径 相关 基因 的 
表达 。 有 研究 表明 叶绿体 内 的 类 胡 莹 卜 素 既 可 以 作为 吸收 光 能 的 辅助 
色素 将 能 量 传递 给 叶绿素 a， 又 可 在 持续 胁迫 条 件 下 将 过 剩 光 能 安全 
耗 散 而 保护 光合 机 构 (Barbara et al., 1996; Liu et al., 2004; Holt et al., 
2005)。 我 们 前 期 的 研究 也 发 现 低温 胁迫 时 甘蔗 能 够 通过 增加 光合 机 
构 过 剩 激发 能 的 耗 散 和 调整 其 叶片 光合 色素 合 量 与 构成 进行 有 效 的 
光 能 利用 和 分 配 。 而 且 本 研究 对 部 分 差异 表达 miRNAs 及 其 靶 基 因 进 


(Т qRT-PCR 1999, AGI 5255825 FZRTUBIRS miR167d 对 靶 基 因 具 
有 负 向 调控 作用 。 因 此 ， 后 续 我 们 对 低温 条 件 下 参与 光合 生理 过 程 的 
miRNAs 进行 深入 研究 其 分 子 调控 机 制 ， 为 抗 寒 育 种 提供 重要 理论 依 
据 和 实验 基础 。 
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